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ABSTRACT 

In itinerant magnetic systems with disorder, the quantum Griffiths phase at T — is unstable to formation of a 
cluster glass (CG) of frozen droplet degrees of freedom. In the absence of the fluctuations associated with these 
degrees of freedom, the transition from the paramagnetic Fermi liquid (PMFL) to the ordered phase proceeds 
via a conventional second-order quantum phase transition. However, when the Griffiths anomalies due to the 
broad distribution of local energy scales are included, the transition is driven first-order via a novel mechanism 
for a fluctuation induced first-order transition. At higher temperatures, thermal effects restore the transition to 
second-order. Implications of the enhanced non-Ohmic dissipation in the CG phase are briefly discussed. 

Keywords: non-Fermi liquid behavior; rare events; Griffiths phases 

1. INTRODUCTION 

When quenched disorder is introduced into a magnetic system in its ordered phase, exponentially rare defect- 
free regions or droplets can form which are nearly magnetically ordered, even if the random magnet itself is 
not. The system is correlated over the length scale of the individual droplets and long time scale (compared 
to that of a single spin) fluctuations are induced due to coherent flipping of these large volumes of spins. Such 
slow fluctuations lead to singularities of the free energy in what is known as the Griffiths phase. 1 Classically, 
these Griffiths singularities are weak, essential singularities 2 and can only make exponentially small corrections 
to thermodynamic quantities. However, near quantum phase transitions, disorder is perfectly correlated in 
imaginary time making its effect much stronger than in the corresponding classical problem. Recently, much 
progress has been made in understanding the effects of disorder on quantum critical points by focusing on the 
dynamics in the quantum Griffiths phase (QGP). 3 ~ 9 

One compelling reason for much the focus in this area is in understanding the role QGP anomalies play in 
the non- Fermi liquid behavior of disordered strongly correlated electron systems. 10 For example, in the random 
transverse field Ising model in dimensions 1, 2 and 3, 5 it has been shown that a quantum Griffiths phase 
characterized by a broad distribution of local energy scales precedes a T — transition in the universality class 
of the infinite-randomness fixed point (IRFP). The distribution of energies follows a power-law form with a 
continuously varying exponent a = d/z' where d is the dimensionality and z' is a dynamical critical exponent 
which diverges at the IRFP. This power-law distribution function leads to strong corrections to Fermi liquid 
scaling in the QGP. In particular, thermodynamic quantities become divergent at T = showing that quantum 
Griffiths anomalies can have significant impact on the behavior of random systems near their critical points. 

Significant progress has been made recently by considering the rare regions in the Griffiths phase as inde- 
pendent droplet degrees of freedom, and an elegant classification scheme was proposed based only on general 
symmetry arguments (for a review see Ref. 11). It was shown that the effects of quantum Griffiths anomalies 
range from weak, classical corrections, to complete destruction (rounding) of the phase transition. However, 
these arguments apply only to insulating magnets for which the droplets can truly be considered as independent. 
By considering interactions induced between droplets in a metallic host, it was recently shown 9 that the QGP is 
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Figure 1. Schematic depiction of the coarse-graining procedure leading to the interacting droplet model of disordered 
magnetic systems. Rare, defect-free regions in the original system (shaded, right hand panel) are identified as the 
important degrees of freedom. Individual droplets then interact via long-ranged interactions mediated by the metallic 
bulk. 



unstable to formation of a cluster glass phase (CGP), and the conventional quantum critical point is destroyed. 
In the present work, we will examine the same model as in Ref. 9 and shown that fluctuations due to Griffiths 
phase anomalies have the further effect to drive the transition first-order. At higher temperatures, thermal effects 
act to suppress these fluctuations and a second-order transition is restored. 

The paper is organized as follows: Sections [2] and [3] introduce the model and formulate the saddle-point 
approximation to it; its solution in the limit of uniform droplet sizes is presented in Section 2] Sections [5j [5] and 
[7] are devoted to solution of the model incorporating the broad distribution of droplet sizes; and Section [8] is for 
summary and discussion. 

2. THE MODEL OF INTERACTING DROPLETS 

By focusing on locally ordered droplets as the relevant degrees of freedom, 7,8 a coarse-grained action in terms 
of a fluctuating, local order parameter can be devised 9 to study the effect of disorder on phase transitions in 
random magnets. This coarse-graining procedure is shown schematically in Figure[TJ In this framework, a single, 
isolated droplet acts as a classical moment with dynamics in imaginary time and can be mapped onto a classical, 
one-dimensional spin chain with a local coupling constant determined by the size of the droplet. The local action 
for an isolated droplet is then given by 

SL,i = f Q drdrV 4 (r)r i (r - r')fa{r') + ^ jf" dr#(r) (1) 

where fa is an iV-component order parameter field for the ith droplet and Ti(r) is a generic coarse-grained 
two-point vertex. For definiteness, we will focus on the itinerant Heisenberg antiferromagnet 12 (N > 1) so that 

Ti(un) = (n + M); (2) 

the w n are bosonic Matsubara frequencies ui n = 2nirT and the local coupling constant is given by r.j w S(Li) d 
where Li is the linear dimension of the droplet, d is the spatial dimension, and 5 < is the coupling constant 
of the clean magnet. Droplet sizes follow a Poisson distribution P(Li) ~ exp[— p(Li) d ] so that, correspondingly, 
P(ri) ~ exp[— pn/S], where p is the volume fraction of droplets. 
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In the absence of interactions between droplets, the action described by |T]) and @ is well understood. 8 
Ohmic dissipation given by the linear uj n term corresponds to long-ranged 1 /r 2 interactions in imaginary time so 
that the corresponding Heisenberg spin chain is exactly at its lower critical dimension 13 and the resulting energy 
gap is found to be 7,8 e, ~ exp[7rr,/«]. This leads to a power-law distribution of local energies scales P(e») ~ e" _1 
where the Griffiths exponent a = d/z' is a non-universal function of parameters and is expected to decrease as 
the magnetically ordered phase is approached. Standard quantum Griffiths phase (QGP) phenomenology follows 
from this distribution; for example, the average susceptibility x ~ T a ~ l diverges as T — > for a < 1. 

The above scenario applies in the absence of interactions between droplets and, in particular, is valid for 
insulating magnets. However, in itinerant systems, electrons in the bulk mediate long-ranged interactions which 
must be included. 9 Formally, the coarse-graining procedure requires that the order parameter fluctuations in 
the metallic bulk be integrated out, leading to an effective RKKY interaction between droplets which has the 
asymptotic form 







S ^ = TJTVa I dr <M T )<M T ) ( 3 ) 
K-tiij) Jo 

where Rij is the distance between droplets i and j, and Jij can be taken to be a random amplitude with zero 
mean and variance (Jfj) = J 2 - The effective description of a random itinerant magnet near a phase transition 
is then given by the full action 

s=E 5 L,i (4) 

i i,j 

Recently, 9 action (|4|) was studied for the Heisenberg antiferromagnet and it was found that infinitesimal 
interactions can destabilize the QGP found in the non-interacting theory, leading to the formation of the cluster 
glass phase (CGP). These results showed that the interactions can have a non-trivial effect on the droplet 
dynamics, even at the saddle-point level which is formally justified in the large- N limit. In the following Section, 
we review the saddle-point theory for the action (g]). 

3. SADDLE POINT THEORY 

By including RKKY interactions between separated droplets, individual droplets no longer have only their own 
dynamics given by the kernel @, but also acquire dynamics through their interactions with an effective bath 
given by the collective fluctuations of all other droplet degrees of freedom. To obtain a description for this 
process, we average over the random variables Jij using the standard method of introducing n copies or replicas 
of the system 14 described by ([4]) and applying the identity 

lynx _ i 

(m^lim^- (5) 

where the partition function Z = Tr{exp[— S]} and (• • -) dis represents averaging over random Jij. The interaction 
term of (j4|) then becomes 

i 72 ™ rP 

%i = -o7^l23 E / drdr'^(r)^(r')^(r)^(r') (6) 

which can be decoupled by introducing a Hubbard-Stratonovich field 15 Qf b (T — r'). The replicated partition 
function can then be written as 

(Z n ) - /D0DQex P {-S eff [0,Q]} (7) 



with effective action 

n 1 n 

s cS i<i>, q] = E E S ^M\ + o E E 



2 

a.b—1 i,LJ-n 



4("n) (8) 



\ E E QfM ((^a) Q?W 



a, 6=1 i,j,LL>n 
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Defining the "cavity field" 

we see that the problem of interacting droplets has been reduced to a new single-site problem, given by 

n 1 n 

SlM^ Ql-EE + 2 E E <^K) A f (io) 

a— 1 i,u? n a,b— 1 i,uj n 

so that the effect of the other degrees of freedom is simply to introduce additional dissipation through the 
frequency dependence of A(u> n ). At the saddle-point level, the cavity field is determined self-consistently from 
the relation 

A ° b K) - E <^-k)^k)> l (ii) 

where the expectation value is taken with respect to Si^cff- In the paramagnetic phase, (0"(w n )^(o; n )) = 
SabXh(fi, ijJn) is diagonal in replica indices. Noting also that the ri are independent of droplet position, is 
equivalent to 

Af(u n ) =gS ab Xh{un) = gS ab / driP(n)xh(ri,Lo n ) (12) 



where 5^ J 2 E ij -(^)" M - 

To finalize the formulation of the saddle-point theory, the correlator Xh can be calculated in the large- N limit 
by decoupling the quartic term in with an auxiliary Hubbard Stratonovich field A^ . The equations governing 
the dynamics of the droplets in the interacting theory then become 
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Xl(^„) = - / dr % ■ — - (13) 



ri + Xi + \u„\- gxU^n) 



Aj = l - _ (14) 

2 u„ r * + X i + - gXL(u n ) 



which must be evaluated self-consistently for Ai and Xh(^n)- 

4. QUANTUM CRITIC ALITY FOR UNIFORM DROPLETS 

To understand the transition in the interacting droplet model, it is instructive to first study equations (fTB"]) and 
(|14|) in the limit of uniform droplet size, P{ri) = 5{r.i — f). In this limit, the problem becomes equivalent to that 
of a quantum spin glass which is known to undergo a conventional second-order transition from a paramagnetic 
Fermi liquid (PMFL) to a spin glass at large enough g. le 

Self-consistency equation (|13[) reduces to an algebraic equation for xi^n) and is easily solved: 



XuK) = (r + A + |w n | ± yj(r + A + \lu h \) 2 - 2gj (15) 

where A = A(f). Clearly, this solution becomes unstable when the energy scale defined by A = (f + A — y/2g) 
becomes zero. The phase boundary can then be determined by setting A = and solving (fT4)) self-consistently 
for r c (T) (or, equivalently, A C (T)). In particular, at T — 0, this gives 

rc(T = 0) = ^ + ^ + ln(#V) (16) 



2tt \2 V 2A , 

where A is an ultra-violet cutoff. In the limit g — > 0, r c (0) — > — oo showing that frustration induced by the 
interactions is necessary to stabilize the glass phase. 
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The phase for A > (f > f c ) is characterized by the frequency dependence of xu- From the solution (fT5|) . 
we find three separate regimes: 

!uj ,lu < A < 

,A<w<V¥ ■ (17) 
xu(o)-^ 1 ,w>V¥ 

The linear, Iow-lj dependence of xu for A > is characteristic of a Fermi liquid, while, right at the transition, 
Xu behaves as y/oJ. This sub-linear form at A = leads to non-Fermi liquid behavior of all thermodynamic 
quantities, though these violations are fairly mild since the corrections remain finite at T = 0. Thus, there is no 
quantum Griffiths behavior in this model since we have neglected the fluctuations arising from a distribution of 
droplet sizes. 

The transition at T = can also be approached from the magnetically ordered side. In this case, the quantum 
critical point can be determined when the mean-field stability criterion for the glass phase vanishes, 17 i.e.: 

A SG = 1- v / 2^Xu(0)=0. (18) 

Using the solution (| 1 5|) for xu, we find Asg ~ vA which vanishes at A = 0. 

Thus, examining the transition through the instability of either the paramagnetic Fermi liquid or the spin- 
glass phase yields a quantum critical point at A = 0, consistent with a conventional second-order transition in 
this problem. In the next Sections, we will consider the model with the full distribution of droplet sizes and 
show that the two instability criteria do not coincide, revealing the singular effect of droplet fluctuations. 



5. DISTRIBUTED DROPLETS AND THE CLUSTER GLASS PHASE 

For dilute impurities, the defect-free regions assume a Poisson distribution which can be written in terms of the 
local coupling constant as 

P( r A = H^ e 27r«(r i -f)/« ; n <f<0. (19) 
u 

The offset f is used to tune through the transition and k — pu/2it6. The uniform limit of the previous Section 
is recovered for k — > oo; however, by accounting for droplets of all sizes, we are also including their fluctuations 
which are associated with Griffiths phase behavior. 



Using the uniform solution (|15|) as a zeroth order approximation, we can iterate equations (|13[) and |T] 
until self-consistency is achieved. Already at one iteration loop (HL), the instability of the Griffiths phase at 
T = is apparent. At this level of approximation, equation (fUT| can be integrated exactly. Identifying the local 
droplet energy = f + + Xi — ffXu(O), we find the relation between energy scale and local coupling strength, 
6j tx exp[2irf~ 1 r i /u] where 

2 ^ A2 +^I (20) 

A + y/2g + \/A 2 + 2Ay/2g 

This allows us to switch from integration over in equation (|13|). to integration over via the replacement 
P(ri)dri — > P(ei)dei where 

P(ci) = (nfe- 27 " if / u ) e nf -\ (21) 



Defining a = nf , (|21[) takes the power-law form characteristic of droplets in the Griffiths phase with Griffiths 
exponent a. The novelty arises when this distribution is used to calculate the frequency dependent cavity field, 

XlK) = f e~ 2 — /" j K de t _f ) — , (22) 

2 Jo e; + M-s(xuM-xu(0)) 

where A € is an upper cutoff in energy. This gives 

X 7 ^-X^M = -lH a - 1 +0(\u;\), (23) 
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(non-Ohmic dissipation) 



Figure 2. Schematic phase diagram of random itinerant magnets. Interactions destabilize the putative quantum Griffiths 
phase at a — 1 , favoring the formation of the cluster glass phase for a < 2. The phase boundary has an exponential tail 
due to the exponentially rare nature of the large frozen droplets. 



where 7 is a constant. Thus, droplets for which a < 2 acquire non-Ohmic dissipation and the corresponding 
Heisenberg spin chains find themselves above their lower critical dimension 13 and can order in the imaginary 
time direction. Due to the frustration caused by the effectively random RKKY interaction between these ordered 
droplets, they form a glassy state dubbed the "cluster glass phase" (CGP). 9 Importantly, the CGP is present 
over the entire range where one would predict the QGP from the isolated droplet model. 8 

The finite temperature phase boundary can also be estimated in this formalism. The critical coupling constant 
is determined from fT4")) at e, = 0: 

u f A 1 

' dcj (24) 

ln(l/7<?) 



tt(2 - a) 

and the number of frozen droplets is rtf r ~ exp(— pri iC /5). The critical temperature is then estimated as 7 

pu 



T c ~ rib ~ exp 



log(l/7<7) 



7r|<y|(2-a) 

which has an exponential tail as a — ► 2. The phase diagram is shown schematically in Figure [2] 



(25) 



6. IRFP OF ISOLATED DROPLETS AT QUANTUM CRITIC ALITY 

It is interesting to note that our 1IL expression for the Griffiths exponent hints at the close analogy between the 
present work and the physics of infinite-randomness fixed points (IRFP). 3 From the expression a = fn with / 
given in equation (|20[) , we would conclude at this level of approximation that the Griffiths exponent vanishes at 
the uniform droplet quantum critical point as 

a - VA. (26) 

This is significant in that it gives a simple access to the IRFP that has been lacking in the past, while capturing the 
physics of the quantum Griffiths phase. Of course, this expression for a is only valid at 1IL, and, in particular, 
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does not hold at the lowest temperatures where the system freezes into the CGP. However, QGP and IRFP 
behavior can still be expected at temperatures above the CGP transition temperature. We will address this 
issue further in the final Section. 



7. FLUCTUATION INDUCED FIRST ORDER TRANSITION 



In the paramagnetic phase with a > 2, the leading behavior of Xl(w) is linear in uj. Writing Xl(^) — Xl(0) 
—m x tu and inserting this into the self-consistency condition (] 14[) gives 



u f A , 1 



ri-f + m(0) = — / dc^— — (27) 

u ( e. L + (1 - m x )A 



2tt(1 - gm x ) 



, 2.(1-^) (2g) 



Analyzing this for r, <C gives 
and, comparing to the 1IL results, defines 

/ = {l-~gm x )- l ] (29) 
a = fn, (30) 

which can be checked against the previous results. Note that expression (|30[) is exact in the paramagnetic phase 
provided that m x = 9xl(^)/(?w| w= o is determined self-consistently, so, following the arguments of Section[5j the 
instability of the PMFL to CGP formation still corresponds to a = 2, to all orders. 

This relation can be further illuminated by calculating the slope exactly from (| : 

_ _1 f (l-~gm x) P(r i ) 

= -2(l-gm x )xl(0). 



Solving this for m x , we get 



l-gm x = (l-2gxt(0))- 1 . (32) 
What is striking about this result is that the term in parentheses on the right hand side of (j32|) is the exact 



analogue of the stability criterion (fT8|) for droplets with distributed site energies; 17 in this case, Asg = 1 — ^9Xl(^)- 
Relating this to the definition (|3"0"|) for a, we can write 

a = X sg k (33) 

so that Asg > at the point where the PMFL becomes unstable (a = 2). This should be contrasted with the 
uniform droplet limit of Section 2] where the instability of the PMFL and the instability of the spin-glass phase 
coincided at a conventional second-order transition. In the present case, the stability criteria signify spinodal 
lines and the CGP and PMFL coexist in a region around a first-order phase transition. 

To gain more insight into the nature of the transition, we recast the problem as an eigenvalue analysis of 
the constitutive free energy for which the solution to (|13p is a minimum. The relevant contribution to the free 
energy can be expressed near the minimum as 

T x = / dtodLo'ix^to) - X oM)r x ( W , W ')(XL(cy) - Xo (uj')) (34) 



so that the eigenvalue determining the stability of the PMFL solution xo satisfies 

x < J dmdu;'(xL(^) -XoH)rx(^^')(XL(t/) -Xo(u')) r „_x 
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Figure 3. Plot of the stability criterion A x versus n _1 for several values of A. The dashed line is the curve at the critical 
A c . Data above this line are for A > A c while data below are for A < A c . The inset demonstrates scaling of the curves, 
as discussed in the text. 



The self-consistency condition (fT3")) for x results from the functional derivative &T x j&x = so that 

dw'T x { U , U ')xL(u;') = xi» - / dr t . (36) 

J n + \i + \w\ - gxh{u) 

and 

' 'dwT x (w J w')xo(« / ) = 0. (37) 
Now, define an eigenvalue A*-™- 1 valid at each step in an iteration scheme 

A( n) < j d^/Q^V) - x^r^XXxiTV) - x£V)) (38) 
x " Jd^xt^-x^^W 



and proceed numerically. From this, the spinodal line defined by X x — can be determined at zero and finite 
temperatures, and compared with the condition Asg = 0. 

The results of this procedure at T — are shown in Figure ([3]). As the transition is approached from the 
PMFL, A x appears to approach a positive value as n — > oo. We identify the critical value of A c at the transition 
as the last point where the iterative method converges; the associated curve A XiC vs. n~ l is indicated by the 
dashed line in the figure and extrapolates to zero at n^ 1 = 0. Below the transition, A x eventually becomes 
negative at some finite iteration step and the procedure breaks down. The inset to Figure ((3]) shows a scaling 
analysis for the same data as in the main panel. Using the ansatz A x /A XjC = A±(|A — Acln*), we find excellent 
scaling with exponent x — 2.35. While the meaning of the exponent is unclear at this point, the scaling shown 
in Figure (O demonstrates that the iterative method is well controlled. 

At finite temperatures, we find that the two stability criteria continue to disagree over a portion of the phase 
boundary away from T — 0, as shown in Figure ([4]). The upper panel shows the boundary line defined by 
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-0.67 -0.66 -0.65 -0.64 -0.63 -0.62 

r 

Figure 4. Upper panel: The phase boundary at finite temperatures is shown by the symbols and corresponds to A x = 0. 
The line is a fit to the data as discussed in the text. Lower panel: The static stability criterion calculated at T c . The line 
is a guide to the eye. Empty squares correspond to first-order transitions, full circles to second-order transitions. 

A x = 0. The solid line is a fit to the data using the form log(T) ~ (f — r c ) _1 ) reflecting the fact that the critical 
temperature is proportional to the number of frozen droplets, as discussed in Section ((5|). The value of Asg at the 
boundary defined by A x = is shown in the lower panel. We can see that Asg > A x = for —0.65 < f < r c (0). 
For f < —0.65, the two criteria merge, indicating that the standard continuous transition is restored. In both 
panels, empty squares represent the first-order phase boundary and full circles represent the second-order phase 
boundary. 

This scenario is consistent with a first-order transition induced by the fluctuations associated with the Grif- 
fiths phase. In the uniform limit, these fluctuations were suppressed since all droplets were of equal size, and the 
transition was of a conventional, second-order type. However, with the full distribution of droplet sizes included, 
the transition at T = was driven first-order. At finite temperature, 1/T acts an effective cutoff in the imaginary 
time direction, suppressing fluctuations of the largest droplets. Thus, the restoration of the second-order transi- 
tion at high enough temperatures is a natural consequence of eliminating fluctuations associated with the longest 
time-scales in the exponential tail of the distribution, and we conclude that these fluctuations are responsible 
for the discontinuous nature of the transition at T — and above. 

Finally, in Figure ([5]), we show the phase boundary on a linear temperature scale and for higher temperatures. 
At low temperatures, there is an exponential tail reflecting the fact that the largest droplets which freeze to form 
the CGP are exponentially rare in this regime, as was shown schematically in Figure 

8. CONCLUSIONS AND DISCUSSION 

The fluctuations of rare, defect-free regions in random magnets can lead to a variety of interesting phenomena near 
quantum critical points. In the present work, we focused on itinerant antiferromagnets with a, T — transition 
between the PMFL and CGP, and saw that the Griffiths anomalies provide a novel mechanism responsible for 
driving the transition first-order. 

This work could be of relevance to experiments on random systems in metallic environments. For example, it 
has recently been shown 18 that singular corrections due to rare regions are necessary in understanding the transi- 
tion in a magnetic field of superconducting thin films in the Little-Parks experiment at T = 0. 19 The non- Fermi 
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Figure 5. The numerically determined phase boundary on a linear temperature scale showing the exponential tail due 
to freezing of rare large droplets. 



liquid behavior arising due to quantum Griffiths anomalies may also be relevant for certain strongly correlated 
materials. Recent studies on such varied systems as the rare-earth intermetallic Tb5Si2Ge2 20 and the colossal 
magnetoresistive Lao.7Cao.3Mn03 21 show power-law diverging susceptibilities reminiscent of quantum Griffiths 
physics. We conclude with a brief discussion of the temperature scales relevant to experimental observations of 
cluster glass physics. 

From equation (|24p . we notice that the effects of non-Ohmic dissipation are only apparent below a crossover 
frequency 

w*=exp[-ln(V75)/(2-a)]. (39) 

Thus, for temperatures below this energy scale, thermodynamic quantities should be dominated by the contri- 
bution due to frozen droplets, giving, for example, x ~ T . Above this temperature scale, Ohmic dissipation 
dominates and the results from the non-interacting droplet model will apply. In particular, thermodynamics 
quantities will diverge as a power-law with exponent a — 1 < as in the quantum Griffiths phase. 
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